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Abstract — We consider a two-way relay network, where two 
source nodes, SI and S2, exchange information through a cluster 
of relay nodes. The relay nodes receive the sum signal from 

51 and S2 in the first time slot. In the second time slot, each 
relay node multiplies its received signal by a complex coefficient 
and retransmits the signal to the two source nodes, which leads 
to a distributed two-way beamforming system. By applying the 
principle of analog network coding, each receiver at SI and 

52 cancels the "self-interference" in the received signal from 
the relay cluster and decodes the message. This paper studies 
the 2-dimensional achievable rate region for such a two-way 
relay network with distributed beamforming. With different 
assumptions of channel reciprocity between the source-relay and 
relay-source channels, the achievable rate region is characterized 
under two setups. First, with reciprocal channels, we investigate 
the achievable rate regions when the relay cluster is subject to a 
sum-power constraint or individual-power constraints. We show 
that the optimal beamforming vectors obtained from solving the 
weighted sum inverse-SNR minimization (WSISMin) problems 
are sufficient to characterize the corresponding achievable rate 
region. Furthermore, we derive the closed form solutions for 
those optimal beamforming vectors and consequently propose 
the partially distributed algorithms to implement the optimal 
beamforming, where each relay node only needs the local 
channel information and one global parameter. Second, with 
the non-reciprocal channels, the achievable rate regions are 
also characterized for both the sum-power constraint case and 
the individual-power constraint case. Although no closed-form 
solutions are available under this setup, we present efficient 
algorithms to compute the optimal beamforming vectors, which 
are attained by solving a sequence of semi-definite programming 
(SDP) problems after semi-definite relaxation (SDR). 

Index Terms — Two-Way Relay, Distributed Beamforming, 
Achievable Rate Region, Pareto Optimal, Semi-definite Program- 
ming (SDP), Semi-definite Relaxation (SDR). 



I. Introduction 

Cooperative communication has been extensively studied 
in past years, where various cooperative relaying schemes 
have been proposed, such as amplify-and-forward (AF) ||T|, 
decode-and-forward (DF) 13, compress-and-forward (CF) ||3l, 
and coded-cooperation H. Among these schemes, due to its 
simplicity, the AF-based relaying is of the most practical 
interest, where multi-antenna relay beamforming has also been 
explored to achieve higher spatial diversity Q. In certain 
resource constrained networks, such as sensor networks, the 
node size is limited such that each node could only mount 
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a single antenna IS). In order to exploit the multi-antenna 
gain in such size-limited cases, distributed relay beamforming 
strategies have been developed where the relaying nodes 
cooperate to generate a beam towards the receiver under sum 
or individual power constraints liUtSJ. 

As an extension to the AF-based one-way relaying scheme, 
the AF-based two-way relaying scheme ||9| is based on the 
principle of analog network coding (ANC) [|10j to support 
communications in two directions. Traditionally, two-way re- 
laying avoids the simultaneous transmissions of two source 
terminals, and requires four time-slots to finish one round 
of information exchange between them. On the contrary, the 
two-way relaying scheme proposed in ifTOl allows the relay 
to mix the data and amplify-and-forward it, where the two 
terminals exploit the underlying self-interference structure. By 
doing so, the amount of required transmission time-slots is 
reduced from four to two and the overall network throughput 
is thus improved. There are several other works discussing 
such two-way relay systems. In particular, the authors in ||9l 
characterized the maximum achievable rate region for the two- 
way relay beamforming scheme by assuming a single relay 
node equipped with multiple antennas and two source nodes 
each equipped with a single antenna. As a counterpart of the 
work in ||9|, the decode-and-forward two-way relaying has 
been studied in ||TT|. The authors in |[T2l studied the AF-based 
two-way relay with distributed beamforming, where the focus 
is to minimize the total transmit power of the source nodes 
and the relay cluster under a given pair of signal-to-noise ratio 
(SNR) constraints. 

The works on characterizing the rate region of two-way 
relaying has also been done in ITSl . where authors considered 
the distributed beamforming case. However, all of existing 
works only obtained numerical solutions. The missing of 
closed-form solutions leads to difficulties in designing efficient 
algorithms due to the lack of insight into the structure of the 
optimal beamforming vectors. Thereby, in this paper, we try to 
seek the closed-form solutions for the optimal beamforming 
vectors to characterize the maximum achievable rate region 
and correspondingly propose efficient distributed algorithms. 
Our work differs from the work in [5:1 from two main aspects. 
First, we assume a cluster of single-antenna relay nodes and 
consider distributed two-way relay beamforming rather than 
the multiple-antenna single-relay beamforming. Due to the 
distributed feature, we will study the case where each relay 
node has an individual power constraint in addition to the case 
where all relay nodes are subject to a sum-power constraint. 
Second, we present closed-form solutions for the optimal 
beamforming vectors when we have reciprocal channels. 
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Fig. 1: System model 



The rest of the paper is organized as follows. In Section HIl 
we introduce the system model and define the achievable rate 
region. In Sections |III] and IIVI we characterize the achievable 
rate regions with two different assumptions on the channel 
reciprocity. Sub-optimal schemes with lower complexity are 
discussed in Section [V] Numerical results are presented in 
Section |Vl] with conclusions in Section IVIII 

Notations: We use uppercase bold letters to denote matrices 
and lowercase bold letters to denote vectors. The conjugate, 
transpose, and Hermitian transpose are denoted by (•)*, (•)^, 
and (•)^, respectively. The phase of a complex variable a 
is denoted as Za. We use tr(-) and rank(-) to represent 
the trace and the rank of a matrix, respectively. A diagonal 
matrix with the elements of vector a as diagonal entries is 
denoted as diag(a). A ^ means A is positive semi-definite, 
a ^ b means a^ > bi component-wise, and stands for the 
Hadamard (elementwise) multiplication. 

II. System Model 

As shown in Fig.lT] we consider a distributed two-way relay 
system consisting of two source nodes S 1 and S2, each with 
a single antenna, and a relay cluster with K single-antenna 
relay nodes i?i's, i = 1, • • • ,K. No direct links between SI 
and S2 exist. The forward channels from SI and S2 to relay 
node i are denoted as /ii,; and /i2.t, respectively, while h\ ^ 
and /12 j denote the backward channels from relay node i to S 1 
and S2, respectively. All the involved channels are assumed to 
take complex values and remain constant during one operation 
period. In addition, all channel state information is revealed to 
SI, S2, and the design/control center where the beamforming 
solution is solved. 

The two-way relaying takes two consecutive equal-length 
time-slots to finish one round of communication between S 1 
and S2 via the relay cluster with perfect synchronization 
assumed among SI, S2, and Ri, i ^ 1,- ■ ■ ,K. In the first 
time-slot, SI and S2 transmit their signals simultaneously to 
the relay cluster; the i-th relay node receives the mixed signal 
ti{n), which is expressed as 



ti{n) = hi^iSi{n) + h2AS2{n) + Vi{n), 



(1) 



where si{n) and S2{n) are the transmitted symbols at time 
index n; and Vi [n) is the receiver noise at relay node i, which 



is assumed to be circularly symmetric complex Gaussian 
(CSCG) with zero mean and variance of. In the second time- 
slot, upon receiving the mixed signal, relay node i multiplies 
a complex coefficient Wi and forwards the signal, which is 
given as Ui{n) — Witi{n). At the source node terminals, SI 
and S2 receive the sum signals from all the relay nodes, which 
are respectively given as 
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hl^^Ui{n) + zi{n), 



/l2.j"i(") + Z2in), 



(2) 



(3) 



where zi{n) and Z2(n) are the noises at SI and S2, re- 
spectively, which are assumed to be CSCG with zero mean 
and variances cr^j and ag^^ respectively. Since SI and S2 
know their own transmitted signals, si(n) and S2{n), re- 
spectively, they could subtract the resulting self-interference 
terms J2i=iK,i'^ihi,i^ii'>T') ™d J2t=i K,i'^^ifT'2AS2{n) from 
the received signals, respectively. Accordingly, the remaining 
signals for SI and S2 are 
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[hl^iWih2aS2{n) + hl^^w^Vtin)] + zi{n){4) 



[h2jWihi,iSi{n) + h2jw^v^{n)] + Z2{n)i5) 



Therefore, for a given w = [wi,--- ,wk]'^ the maximum 
achievable rates for the end-to-end link from S2 to SI and 
from SI to S2 are respectively given as 
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hi h^, f2 
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= h2 h^ with hi 
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(7) 
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[hi^i,--- JkkV and h[ = [Kir--,hlj^] 
In addition, Ai = diag[\hl^^\'^af,--- ,\K,k\'^'^kI ^2 = 
diag[|ft,^_ipo-2 ... ,\h^j^\'^a%], Psi and P52 are the maxi- 
mum transmit powers at S 1 and S2, respectively, and the factor 
1/2 is due to the use of two orthogonal time-slots for relaying. 
Accordingly, we can define the set of rate pairs achievable 
by all feasible beamforming vector w's as 



7^= U {{ri,r2):ri <Ri,r2<R2}, 



(8) 



weOu 



where the feasible set fiu, can be defined by either a sum- 
power constraint or individual-power constraints. Specifically, 
when the sum-power constraint is considered, we have fluj ~ 
{w : pr{'w) < Pr}, where Pr is a scalar power limit 
and pr{w) is the sum-power of the relay cluster given the 
beamforming vector w. When individual-power constraints are 
considered, we have il^ = {w : P7j(w) ^ P_r}, where 
Pfl(w) is a vector of individual transmit powers, Pr is a 
vector with its elements denoting the power constraints for 
individual relay nodes, and ^ is element-wise. 

When time-sharing between different achievable rate pairs 



is considered, the achievable rate region is then defined as the 
convex hull over the set of TZ. 

Definition 1: The achievable rate region O is the convex 
hull over the set of achievable rate pairs TZ, i.e., 



= i7evx(7^), 



(9) 



where i7cvx(-) is the convex hull operation. 

The goal of this paper is to efficiently characterize the 
achievable rate region O. According to different assumptions 
on the channel reciprocity between the forward and backward 
channels, we first study the reciprocal case, and then study the 
non-reciprocal case. 

III. Reciprocal Channel Case 

In this section, we assume that the forward channels from 
each source node to the relay nodes are reciprocal to the 
backward channels from the relay nodes to each correspond- 
ing source node, i.e., /ii ^ = K^ ^ and /i2,i = K^ j, for 
i = 1,- ■ ■ ,K, which usually holds for a time-division-duplex 
(TDD) relaying system. In this case, it is obvious that when 
Zw = — (Zhi + z^h2), both rates given by (|6]l and (|7J are 
maximized for a given set of jwil's. Thus, we only need to 
further find the optimal amplitudes for the elements in w. 
Let Xi = \wi\ and fi = |/ii,i||/i2,i|; we rewrite (|6]l and ^, 
respectively, as 



Ri 



R2 




(10) 



(11) 



Av^kW 



and 



1 |^i,k||^2,a'|]^- In order to obtain O, 



where x = [\wi\ 

[\hi,i\\h2.i\,--- , \hi,K\\h2,K 

we need to characterize the Pareto boundary of TZ. A 

common method is via solving a sequence of weighted 

sum-rate maximization (WSRMax) problems lfT4l . each for a 

different non-negative weight vector (A, 1 — A), < A < 1, as 

follows 
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(12) 
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Unfortunately, we cannot derive the closed-form solution for 
the WSRMax problem. However, from ( fTOl i and (fTTT i. we see 
that the received SNRs at SI and S2 are 

Ps2{f^? 



SNRi = 



SNRa = 



(t|, + x^Aix' 



Psi{{ x) 



(14) 



(15) 



"S2 + X" A2X 

respectively, where their numerators differ by only a scalar 
constant. As shown later, for each given weight vector 



{fi,l — IJ.),0 < fj, < 1, we could easily find a closed- 
form solution for the following weighted sum inverse-SNRs 
minimization(WSISMin) problem 



'SI 



s.t. 



x^Aix 

^ Ps2(f^x)2 



+ (i-m)' 
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Psi(f^x)2 



(16) 
(17) 



Hence, we could quantify the Pareto boundary for the inverse- 
SNR region. Based on this observation, together with the fact 
that there exists a bijective mapping between an inverse-SNR 
pair and a rate pair, we are inspired to probe the question on 
whether we could construct the achievable rate region O from 
the easily obtainable inverse-SNR region. In the following, 
we first introduce some definitions related to the inverse- 
SNR region and then show that we indeed can construct the 
achievable rate region O from the inverse-SNR region, based 
on the set of closed-form solutions for a sequence of the 
WSISMin problems. 

A. Characterizing the Achievable Rate Region 

At first, we introduce some definitions. 

Definition 2: Consider a bijective mapping U : (x,y) 1— >• 
(i log2(l + 1/x), \ log2(l + 1/y)) with (x, y) G R\+; then 
the set of achievable inverse-SNR pairs X is defined as 

X = {{t^M) ■U{tiM) ^Tl} ■ 

For regions TZ and I, we are particularly interested in their 
Pareto boundaries, which are defined as follows. 

Definition 3: The Pareto boundary of 7?^ is defined as 
7' = {(rl,r2):(rl,r2)e7^,((rl,r2) + K)n7^ = (rl,r2)}, 
and the Pareto boundary of I is defined as B = {(^1,^2) : 
ih,t2) eI,iih,t2)-K)f]X = (ii,t2)}, where K = R^ 
is a non-negative cone 1141 . 

Definition 4: Define the points obtained by solving the 
WSISMin problem with a given weight vector (/i, jT) as 

Si^l,X)^{(tl,t2)■. min ^ltl+fit2}. (18) 

(ti,t2)ei 

where p. = 1 — ji. 

Definition 5: The set of points that can be obtained from a 
sequence of WSISMin problems is given as 



0</j<l 



(19) 



In order to show that we can construct O from S{X), where 
O can be obtained by convex hulling over T', we need to prove 
two things: (1) A point in B could be mapped to a point in 
T' and vice versa, which means U{B) = T^; (2) The points 
in B that cannot be obtained by WSISMin are mapped to the 
points in T^ that are unnecessary for constructing O by convex 
hulling over T', i.e., U{B\S{X)) C T'non, where T'non denotes 
the points in T' that are unnecessary for constructing O with 
convex hulling over T'. With the above two statements hold, 
it is easy to see that T' \ T'non C U{S{X)), i.e., the points 
obtained by WSISMin suffice to construct O. 

Proposition 1: The Pareto boundary of the inverse-SNR 
region can be mapped to the Pareto boundary of the rate region 



7?^ by mapping U as given in Definition |2] and vice versa, i.e., where 



M{B). 



(20) 



Proof: See Appendix in IVIII-AI ■ 

Proposition 2: The image of S \ S{I) is not necessary for 
constructing the achievable rate region O. 

Proof: See Appendix in |VIII-B| ■ 

Theorem 1: The points in S{I) are sufficient to construct 
the achievable rate region O. 

Proof: Since V = U{B) and U{B \ S{I)) in V is not 
necessary for constructing the achievable rate region O, it is 
easy to see that U{S{I)) suffices to construct the achievable 
rate region O given that O is obtained by convex hulling over 
V. ■ 

Since we have shown that S{I) suffices to construct the 
achievable rate region O, instead of studying the problem in 
( fT2] i and ( fT3T l, we now study the solutions of the WSISMin 
problems in the following. 



B. Distributed Beamforming under Sum-power Constraint 

In this subsection, we consider the case where the relay 
cluster has a sum-power constraint. The total transmit power 
of the relay cluster is 



PR 



K 

E 

x^Dx 



{\x^h^A^Psi + \x,h2A^Ps2 + |x,|V2)(21) 



(22) 



where we have D = diag[|/ii^ipP5i + |^2,iP^S2 + 
cr?,--- Ahi,K?Psi + \h2,K?Ps2 + (jV\- According to the 
discussion in the last subsection, to quantify the rate region is 
equivalent to seeking the optimal solutions for the WSISMin 
problems given the sum-power constraint as follows: 



min /^/SNRi + /l/SNRs 

X 

s.t. x^Dx < Pr, 



(23) 
(24) 



where < /i < 1 is the weight. First of all, the optimal 
X* must satisfy x*^Dx* = Pr; otherwise, we can always 
scale up x such that the objective function is decreased. 
When /I = or 1, this problem degrades to find the optimal 
beamforming vector for distributed one-way relay, which has 
been extensively studied |5]|7J|151. 

Given the SNRs from ( fT4l i and ( fTSl l and the fact of 
x*^Dx* = Pr, the optimal x for (|23] | should be the solution 
of the following problem: 

. x^[^D/P;,^ + m/^siAi + m/Ps2A2]x 

mm 7^ , (25) 

x^frx 

where v = A'cr|]^/Ps2 + J^<^'s2l-Psi- The above problem is 
equivalent to 



x^fFx 



max ■ 



X X^[l.D/P;^ + /i/^5lAi + /2/Ps2A2]x' 

where the optimal solution is given as 

x* = cr-if/||r-ifi|, 



(26) 



(27) 



diag 


Pi ^ Pk ^ 


, (28) 


(tI + Psi\hiA'' + Ps2\h2A^ , 


(29) 


-H\ 


/^l,»|V/^Sl + |/i2..|'/i/^S2 


) , (30) 



'7j = ^ 

and ^ is a scalar such that x*-^Dx* = Pr. By searching over 
all /i's, we derived a set of x*'s and hence we could compute 
a set of rate pairs by injecting jT/j i into (fTOl i and ( fTTT i. The 
achievable rate region O is then obtained by convex-hulling 
over such a set of rate pairs. 

Partially distributed implementation: The control center 
first decides the appropriate /i such that SI and S2 achieve 
a desirable rate pair; and it broadcasts /i and the global 
constant ^/||r^^f||, while Psi, Ps2, fsi, 0-32 are constant 
and assumed to be known at all the relays. Upon receiving 
the broadcast message from the control center, each relay node 
determines the optimal Wi from its local information hi,i and 
/i2,i, which is given as 



w. 



e 



\h1A\h2 



|r-if|i i^Pi/PR + vi 



-jUhu,+Zh2,0 



(31) 



C. Distributed Beamforming under Individual-Power Con- 
straints 

In the previous subsection, we assume that the relay cluster 
has a sum-power constraint. In practice, each relay may have 
its own power constraint due to the individual power supplies. 
The transmit power at relay i is given as 



PRS 



\\hiA^Psi + \h2A^Ps2+a^), (32) 



where p^ ,; < pi, with pi is the maximum allowable power 
for relay node i. Equivalently, we could set p^i = afpi with 
< c^i < 1 as a new design variable. Correspondingly, the 
received SNRs can be rewritten as (|33] | and (l34l i 
Let 



diag[alpi\h 



1,11 



diag[cr^pi|/i2,iP, 



Hi 
H2 

gi = \/Pt\hl,i\\h2,i\/y^Di,i. 

We can recast (l33T l and ( |34] | as 

Ps2a^Sg^a 



,(^]iPK\h2,KW 



SNRi 
SNR2 



fT|2 + a^H2D-ia' 



(35) 
(36) 
(37) 



(38) 
(39) 



respectively, where ^ a ^ 1. The WSISMin problem for 
the individual-power constraint case is now given as: 

^ a^(HiD- V/^S2 + H2D-i/2/P5i)a 



mm ■ 
o-<Q-<i 



a'^gg^a 



which is equivalent to solve 
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+Psi\hi,i\^+Ps2\h2,: 
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For notation simplicity, let 

* = [(HiA*/Ps2+H2A/Psi)D-Vi.]'/' 



g 



g 



(42) 
(43) 



where ^ is diagonal with its diagonal elements denoted as ^i, 
i = 1, • • • ,K. Then the above problem becomes 

\2 



{S,a) 



O^a^l 1 



*a 



(44) 



For each given fi, (l44l i can be solved analytically by following 
the results in m. Before we present the solution, we first 
define (j>i = Qi/ipf for i = 1, • • • , A' and (/)a'+i = 0- Then we 
> ■■■ > 0T/<- > 4>TK+i- Moreover, let 



'i as (j)r^ > 



sort ( 

a*-*^^ as 

^ \ Afe0j, j = Tfe+l,-- -TK 

Then the solution for ( |44] | is given by following theorem. 

Theorem 2: The solution of (|44] | is a''^ ' given by ( |45] l, 
where fc* is the smallest fc such that A^ < 



and define the j-th element of the vector 



(45) 



given 

Tfc+r 



Proof: This result directly follows the results in fH. ■ 

Partially distributed implementation: Besides the value of 
/i, the control center only needs to broadcast Afe» at each 
operation period. Each relay node then determines 0^ with 
its local information. If 07^ — ^fc*' ^^ relay node transmits 
at its maximum power. Otherwise, it transmits with power 



ct] 



-j(Z/li.i+Z/l2,,) 



{Xk'(j>i)'^Pi, i.e., the optimal w. 

where a 

in general some relay nodes may not transmit with maximum 

transmit power. 



/Pie 
is given in (1451 ). From the solutions, we see that 



IV. Non-reciprocal Channel Case 



In the last section, we have discussed the case where the 
uplink and downlink channels are reciprocal. In this section, 
we discuss the case where the uplink and downlink channels 
are non-reciprocal, which may be the result of deploying 
frequency-division-duplex (FDD) system. 

Due to the lack of channel reciprocity, the approach taken in 
the last section does not apply here. In order to characterize the 
boundary of the region TZ, as we discussed before a commonly 



used method is to solve the following 
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W^A2W 



(46) 
(47) 



for each given weight vector (A, 1 — A). However, the above 
problem is non-convex since the objective function is not 
a concave function. To efficiently quantify the rate region, 
here we resort to an alternative method called the rate-profile 
method |f9l, formulated as 



max 



S.t. 
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- log. 



w e rij, 
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w^Aiw 
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> KRsumi49) 

> RRsurni50) 

(51) 



where Rsum is the sum rate given a rate profile vector [k, k] 
with < K < 1 and K = 1 - K. Let Fi := tTtT, F2 



f*fT 

^2 '■2 ' 



and X = WW . The above problem is equivalent to 
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Pg2tr(F2X) 

(7|i + tr(AiX) 

Psiti-(FiX) 
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(52) 
U53) 



42 + tr(A2X) 



x^o, 

rank(X) = 1, 



> fiRsurni54) 

(55) 
(56) 
(57) 



where the last constraint rank(X) = 1 comes from the fact 
X = ww-^, and Ox = {X : X = ww^jW e ilw} and fi^ 
is defined after dS). According to different assumptions on the 
power constraint, the above problem can be further converted 
into different semi-definite programming (SDP) problems after 
semi-definite relaxation (SDR). 



A. Sum-power Constrained Case 

In this subsection, we assume that the relay cluster operates 
under a sum-power constraint Pr. Given the sum-power 
constraint, the power constraint in ( fSSl l can be replaced by 
tr(DX) < Pr, where D = diag[|/ii,i|2Psi + |/i24pi^S2 + 



^i, • • • , |/ii,KpPsi + \h2,K?Ps2 + 0-i-]. Since the rank-one 
constraint is not convex, the problem is still not a convex 
problem and hence may not be efficiently solvable. To address 
this issue, let us first remove the rank-one constraint and 
consider the following relay power minimization problem for 
given set of k and Rsum = r: 



mm 

X 



s.t. 



tr(DX) 

^ log2 ( 1 



2 
Xh 0, 

which is equivalent to 



Ps2tr(F2X) 
cr|^ +tr(AiX) 

Psitr(FiX) 
42 + tr(A2X) 



> Kr, 



(58) 
(59) 



mm 

X 



s.t. 



tr(DX) 

Ps2tr(F2X) 
cr|i +tr(AiX) 

Psitr(FiX) 
ct|2 + tr(A2X) 
Xt 0, 



> 



> 



71, 



72, 



> Sr, (60) 
(61) 

(62) 
(63) 

(64) 
(65) 



mm 

X 



s.t. 



where 71 = 2^'^''— 1, 72 = 2^'^'' — 1, and they can be considered 
as the SNR constraints for SI and S2, respectively. Since cr|j^ + 
tr(AiX) > and ag2 + tr(A2X) > 0, we could rewrite the 
above problem as following SDP problem: 

tr(DX) (66) 

tr[(Ps2F2 - 7iAi)X] > 71 a^, (67) 

tr[(PsiFi - 72A2)X] > 72a2, ^^g^ 

X ^ 0. (69) 

Denote the optimal value of the above problem as p|j, which 
is the minimum sum-power required by the relay cluster to 
support the target SNRs 71 and 72 for S 1 and S2, respectively. 
If P*R < Pr^ then (71,72) must be an achievable SNR 
pair. Otherwise, 71 and 72 are not achievable. Based on this 
observation, we propose the following bi-section algorithm 
such that the problem ( |52] | without rank-one constraint can 
be solved by solving a sequence of convex power feasibility 
problems, with the assumption that we know an upper bound 
for Rsum, denoted as r^nax- 
Algorithm 1: 

• InitiaUze riow = 0, r„p = rmax- 

• Repeat 

1) Set r ^ ^{riow +rup)- 

2) Solve problem (|66]|-(|69]| with the given r. 

3) Update r with the bi-section method HH: If pjj < 
Pr, set rioui ~ r; otherwise, r„p = r. 

• Until Tup — riow < e, where e is a small positive accuracy 
parameter 

The rate upper bound r„iax can be derived as follows. 
We first decouple the two-way relay channel into two one- 
way relay channels and obtain a rate for each one-way relay 
channel. Denote the larger rate as f. Then r^ax can be set 
as 2f. The one-way distributed relay beamforming with sum- 



power constraint is well-studied, and the rate can be derived 
from the results in Q. 

1) Rank-one solution: The resulting optimal solution 'Kopt 
obtained from Algorithm 1 may not be of rank-one due to 
the SDP relaxation, which means that Xopt may not lead to 
an optimal beamforming vector w. However, since there are 
only two linear constraints (|67] | and (l68l l. it has been shown 
in 191 and lfT6l that an exact rank-one optimal solution can 
always be constructed from a non-rank-one optimal solution. 
The transformation techniques developed in |9| and |fT6l 
can be used to obtain the rank-one solution. Note that the 
beamforming solution for the non-reciprocal channel case is 
fully centralized, which cannot be implemented in a partially 
distributed fashion. 



B. Individual-Power Constrained Case 

In the previous subsection, we have discussed the sum- 
power constrained case where the non-convex rate maximiza- 
tion problem is converted into a sequence of convex sum- 
power minimization problems. In this subsection, we put a 
stricter limitation on the relay power by assuming that each 
node has its individual power constraint. In this case, following 
a similar SDR technique to that in the previous subsection, the 
optimization problem with individual power constraints can be 
cast as 



max 

X,_Rs„™ 



S.t. 



Pg2tr(F2X) 
cr|i +tr(AiX) 

PsMFiX) 
a|2 + tr(A2X) 

X^O, 



> 7i, 



> 



72, 



,R\ 



(70) 
(71) 

(72) 

(73) 
(74) 



where D,;.i and X^.; are the i-th diagonal elements of D 
and X, respectively. The transmit power at node i amounts 
to Di iXi i and the individual power limit at node i is 
PJj.. However, we cannot translate the above problem into a 
sequence of power feasibility problems as given in the last 
subsection, since we now have K individual power constraints 
rather than a single sum-power constraint for the whole relay 
cluster. Alternatively, we aim at solving a sequence of the 
following problem via bi-section search over r. 



max 

X,r 



S.t. 



r (75) 

tr[(P52F2 - 7iAi)X] > 7ia?, (76) 

tr[(PsiFi - 72A2)X] > 72^2, (77) 

X{i,i)<PyD{i,i), i = l,--- ,K, (78) 

X ^ 0. (79) 

The above problem is convex over X at each given value of 
r. Let r* be the maximum value obtained by solving ( iTSl l. For 



a given value of r, we solve the following feasibility problem that has the minimum v{wi), i.e.. 



Find X (80) 

s.t. tr[(Ps2F2 - 7iAi)X] > jial (81) 

tr[(PsiFi - 72A2)X] > ^2(jI (82) 

X(z,z)<PyD(z,z), ^ = 1,•••,A^ (83) 

X ^ 0. (84) 

If it is feasible, we have r < r* and the corresponding rate is 
achievable. Otherwise, we have r > r* and the corresponding 
rate is not achievable. Based on this observation, we apply 
bi-section search over r to solve the problem in (fTsT l, where 
we solve a convex feasibility problem of (|80] l at each step. We 
start with an interval [0, r,nax\ that contains the optimal value 
r* where rmax can be obtained in a similar way as that for the 
sum-power constrained case, and run the following algorithm. 
Algorithm 2 

• InitiaUze riow=0, r„p = r^ax- 

• Repeat 

1) Set r ^ ^(rioni + ?>). 

2) Solve the feasibility problem given by (l80ll-(l83]l 
with given r. 

3) Update r: If the problem is feasible, set rio^ = r; 



otherwise, r. 



Until r. 



up 



rio 



up — ' ■ 

< e. Then r* 



rio 



1) Rank-one solution based on randomization: Similar to 
the sum-power constrained case, the solution of X at the 
end of Algorithm 2, denoted as 'S.opt, may not be rank- 
one. However, since there are K + 2 linear constraints here, 
we cannot apply the rank-one decomposition technique in 
lfT6l . which require the number of linear constraints to be 
less than or equal to 3. Fortunately, various techniques have 
been developed ifTTl to generate good rank-one approximate 
solutions to the original problem. One such efficient approach 
is based on randomization iflTl : using Xopf to randomly 
generate a set of candidate weight vectors, {w;}, from which 
the "best" solution for the beamforming vector w is selected. 
There are three ways of generating {w;} as presented in 
ifTTl . In order to satisfy the individual power constraint, we 
adopt the routine named randB in iflTl . Specially, let e; be 
the vector whose elements are independent random variables 
uniformly distributed on the unit circle in the complex plane, 
i.e., its i-th element [e;]i — e^^'S where 6';.i's are independent 
and uniformly distributed over [0, 27r). We choose w; such 
that its i-th element [w;]i = ■^/[XoptJi.Je;]^. As we see, 
|[w;]ip = \JLopt\ii', hence the individual power constraint can 
be satisfied. 

For each X'^')=Wiw/^, we associate each \^i with a value 

w(W(), 



v{^i) 



PsiFi 



72^2 



max I 1 — tr[( 

1-tr 2 — 

71 ^Si - 



a 



-)w,wf] 



S2 



1 N H- 

2-)W(W; 

SI 



(85) 



which reflects how much the constraints are violated. The 
"best" weight vector among the candidate vectors is the one 



/* 



w = 



argminw(w;), 



(86) 
(87) 



V. Sub-optimal Schemes 



In this section, we propose some suboptimal schemes with 
lower complexity for implementation than the optimal ones 
established in the previous sections. 

A. Reciprocal Channel Case 

In the reciprocal channel case, at first the transmit phases 
OiS, at the relays are matched to the channels as Oi = 
— (Z/ii,i + Zh2.i)- Then with the sum-power constraint, we 
propose the sub-optimal equal power beamforming where each 
relay transmits with equal power. With the individual-power 
constraints, we propose the max-power beamforming where 
each relay transmits with its maximum power. 

1) Equal-power beamforming: All the K relay nodes trans- 
mit with the same power Pfj/K; 6'i's and x/s for 
i = 1, • • • , /v, are given as: 



= -(Z/li,i + Zft.2,j), 



(88) 
(89) 



\ K{Psi\h,,^\^ + Ps2\h2A^ + af) 

2) Max-power beamforming: Each relay transmits with its 
maximum allowable power PR^f, 9i's and x^'s for i ~ 
1, • • • ,K, are given as: 



= -(Z/li,i + Z/l2,j), 



Pr,^ 



Psi\hi.,\^ + Ps2\h2.\^ + a' 



2 I ^2- 



(90) 



(91) 



These sub-optimal schemes enjoy implementation simplicity 
since each relay only requires the local channel information 
hii and h2.i to decide the transmit phase and Xi. 

B. Non-reciprocal Channel Case 

For the non-reciprocal channel case, since the transmit 
phase cannot be matched to the two-directional channels 
simultaneously, we propose a sub-optimal scheme that greedily 
chooses the transmit phases. Specifically, each relay chooses 
the transmit phase to be either Z/ii.i + ZK^ j or Z/i2.i + ZLh\ j, 
whichever maximizes its own contribution to the overall SNRs 
at SI or S2 without considering any other relays' contributions, 
i.e., the transmit phase for each relay is chosen by following 
criterion: 

^ xjPs2\h2.hl,fi^'^\' x2Fsi|/ii,«/iL.e^'f \ 



= argmax 



^l, + ^hK,V^^ ' ^l2 + ^\K^^ 



(92) 

where Xi is the transmit amplitude. To determine x^'s, we pro- 
pose equal-power beamforming for the sum-power constraint 
case and max-power beamforming for the individual-power 
constraint case, which are given in (|89] l and ( |9T| l, respectively. 
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Fig. 2: Achievable rate regions for the reciprocal channel case, 
Psi=Ps2=0 dB, PR=iO dB. 



Fig. 4: Achievable rate regions for the sum-power constraint case, 
Psi=Ps2=0 dB, Pr=10 dB. 
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Fig. 3: Achievable rate regions for the non-reciprocal channel case, 
Psi=Ps2=0 dB, Pfl= [2.5, 3, 0.5, 1, 3] W. 



Fig. 5: Achievable rate regions for the individual-power constraint 
case, Pgi=Ps2=0 dB, Pr= [2.5, 3, 0.5, 1, 3] W. 



VI. Numerical Results 

In the section, we present numerical results to quantify the 
achievable rate region for the two-way relay network with 
distributed beamforming. We assume that the relay cluster 
consists of 5 nodes; the channel coefficients /ii ^ and h\ ^, i ~ 
1, • • • ,K, are independent CSCG variables with distribution 
CAf{0, 1); the channel coefficients /i2,i and /ij j, i = 1, • • • ,K, 
are also independent and distributed as CM{0, 1). The noises at 
the relays and source nodes are assumed to have unit variance 
in the simulations. We change /i from to 1 with step 0.1 and 
obtain 11 Pareto boundary points. For each point, we run 100 
channel realizations to measure the expected performance. We 
then do convex hulling over these points. 

First, we investigate the achievable rate region when chan- 
nels are reciprocal, where we set hi^i = h^ ^ and /i2,i = ^2 i 
for i = 1, • • • ,K. Fig.|2]shows the achievable rate regions with 
the sum-power constraint and individual-power constraints, 
respectively. For the sum-power constraint case, the relay 
power Pji — 10 dB (dB is relative to the unit noise power) 
while the transmit powers Psi = Ps2 = dB. For the 
individual-power constraint case, the relay power constraints 



are given as 2.5, 3, 0.5, 1, 3 W, which is summed up to 10 
dB. For the non-reciprocal channel case. Fig. [3] shows the 
achievable rate regions with the sum-power and individual- 
power constraints, respectively. The powers are the same as the 
settings in the reciprocal channel case. We use CVX, a Matlab- 
based optimization software ITSll . to solve the SDP problems. 
As we see in both Fig. |2] and Fig. [3] due to the symmetry of 
the transmit powers and channel statistics, the achievable rate 
region O is symmetric. When Ps2 = 0, the rate pairs collapse 
to the segment on the horizontal axis, which corresponds to 
the achievable rate for a one-way relay network where only 
SI transmits. Moreover, the rate region for the individual- 
power constraint case is smaller than that for the sum-power 
constraint case. This is quite intuitive since the individual- 
power constraint is stricter than the sum-power constraint. 

In Fig. |4] and Fig. [5j we compare the rate regions for the 
reciprocal and non-reciprocal channel cases under the same 
power constraint assumption, where Fig. |4] is for the sum- 
power constraint case and Fig. |5] is for the individual-power 
constraint case. In both Fig. |4] and Fig. |5] we can see that 
the maximum rate for S 1 in the reciprocal channel case is the 
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Fig. 6: Achievable rate regions and equal-power beamforming rates 

with a sum-power constraint, reciprocal channel, Psi=Ps2=0 dB, 

Pr=0, 10, 20 dB. 



Fig. 7: Achievable rate regions and max-power beamforming rates 

with individual-power constraints, reciprocal channel, Pgi=Ps2=0 

dB, total power=0, 10, 20 dB. 



same as the one in the non-reciprocal channel case. This is 
because such a maximum rate is obtained by optimizing the 
one-way link from SI to S2 without considering the link from 
S2 to SI. Since the one-way link from SI to S2 consists of hi 
and hj, whether h^ = hi or not does not affect the statistics 
of the one-way link from SI to S2. The same argument holds 
for maximum rate at S2. We also observe that the rate region 
for the reciprocal channel case is larger than that in the non- 
reciprocal channel case given the same settings of powers and 
noises. The reason is that we can match the beamforming 
phase to the overall channel phase (i.e., Zw — Zhi + Zh2.) in 
the reciprocal channel case, while we are not able to do so in 
the non-reciprocal channel case. Therefore, TDD based system 
is more favorable in terms of the achievable rate region if the 
channel coherence time is larger than one operation period 
and the transmit-receive chain calibration f 19] can be properly 
done. Besides the rate region, the amount of information 
needs to be broadcast by the control center is significantly 
different. In the reciprocal channel case, the control center 
only needs to broadcast one scalar at each time slot. However, 
in the non-reciprocal channel case, the control center needs to 
broadcast the beamforming vector, which is a complex vector 
of dimension K. 

At last, we investigate the performance of the sub-optimal 
schemes in relative to the maximum achievable rate regions. 
As we see in Fig. |6] for the sum-power constraint case with 
reciprocal channels, the rate pairs achieved by the equal-power 
beamforming scheme, denoted as single points are strictly 
sub-optimal. On the contrary, as shown in Fig. |7] for the 
individual-power constraint case with reciprocal channels, the 
rate pair achieved by max-power beamforming gets closer to 
the boundary when the power budget is reduced!^ In Fig. 
m and Fig. |9] we consider the non-reciprocal channels and 
show the performance of equal-power beamforming and max- 
power beamforming with greedy phase selection as given 



'We set the individual powers P^ = [2.5, 3, 0.5, 1, 3] W with total power 
equal to 10 dB. When total power is changed to and 20 dB, we scale the 
vector proportionally. 
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VII. Conclusion 

In this paper, we considered the two-way relay networks 
with distributed beamforming and investigated the achiev- 
able rate region, which is defined as the convex hull of 
all achievable rate pairs. We studied both the reciprocal 
and non-reciprocal channel cases. In the reciprocal channel 
case, we characterized the rate region when the relay cluster 
is subject to either a sum-power constraint or individual- 
power constraints, respectively. It was shown that we could 
characterize the whole achievable rate region via the Pareto- 
optimal beamforming vectors obtained from solving a se- 
quence of WSISMin problems. Furthermore, we derived the 
closed-form solutions for those optimal beamforming vectors 
and consequently proposed partially distributed algorithms to 
implement the optimal beamforming, where each relay node 
only needs its own local channel information and one global 
scalar sent from the control center. For the non-reciprocal 
channel case, we used the rate-profile approach to compute the 
Pareto-optimal beamforming vectors. When the relay cluster is 
subject to a sum-power constraint, we computed the optimal 
beamforming vector via solving a sequence of relaxed SDP 
power minimization problems followed by a special rank-one 
reconstruction. When the relay cluster is subject to individual- 
power constraints, we solved a sequence of relaxed SDP 
feasibility problems and the rank-one solution is obtained 
by randomization techniques. From the numerical results, we 
found that the achievable rate region is larger in the reciprocal 
channel case than that in the non-reciprocal channel case. 
Hence, TDD-based relaying scheme is more favorable for the 
two-way relay network with distributed beamforming. 



VIII. Appendices 



Fig. 11: Achievable rate regions and equal-power beamforming 

rates with a sum-power constraint, non-reciprocal channel, 

Psi=Ps2=0 dB, Pr=0, 10, 20 dB. 



in ( |92| l. The performance of both equal-power beamform- 
ing and max-power beamforming schemes degrades as Pr 
increases. Thereby, the sub-optimal schemes for the non- 
reciprocal channel case works well only when Pr is small. 
In addition, we also evaluate the performance of equal-power 
beamforming when channel statistics are asymmetrical, where 
we set hi^^,h\i - CNiO.l) and /J2,i,/i2,i -^ C7V(0,i) for 
i = 1, • • • ,K. For the reciprocal channel case, we further set 
h2.i = ft^2 i ^^^ ^'i,i ~ ^2 i- ^^ '^his casc, the equal-power 
beamforming scheme as shown in Fig. [TO] and Fig. [TT] has 
a larger gap to the boundary than that in the case where 
channel statistics are symmetric (i.e., /ii i, h\ ,, /i2.i and K^ i ^ 
CJ\f{0, 1), for i = 1, • • • , A'.), which are shown in Fig. |6] and 
Fig. [8] for the reciprocal and non-reciprocal channel cases, 
respectively. This is due to the fact that the average transmit 
powers for different relays with optimal beamforming should 
appear non-uniform if the channel statistics are asymmetric. 
Therefore, the gain attained by optimal beamforming becomes 
more significant in this case. 



A. Proof of Proposition Q] 

Proof: We will show this by contradiction. Assume 
(a, 6) e B but U{{a,b)) G V. Then we can find another 
point {c,d) E TZ such that c > l/21og2(l + 1/a) and 
d > l/21og2(l + 1/b). According to the definition of I, 
the point ( gs^irf , 22731 ) ^ ^- Thus, there exists a point in 
X such that ^^h^ < o, and js^^rr < ^' which contradicts 
the assumption that (a, b) is a Pareto optimal point. Hence 
U{B) C v. The converse that U{B) 3 V can also be proven 
in the similar way. Therefore, V = U{B). ■ 



B. Proof of Proposition |2] 

In order to prove Proposition |2l we first introduce the 
following two lemmas. 

Lemma 1: Suppose q{x) is a positive, decreasing, and linear 
function with x > 0. The bijective mapping U maps (x, q{x)) 
to {y,p{y))', then p{y) is a non-negative, decreasing, and 
convex function. 

Proof: Let y = log2(l + 1/a;) and p{y) = log2(l + 
l/q{x)) be an implicit function of y, where x > 0. Since 
q{x) is positive, decreasing, and linear, we have q{x) > 0, 
q'{x) < 0, q (x) ~ 0, and hence p{y) > 0. The first-order 
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dp'{y)/dx 
dy/dx 



In 2 



q" {x){x + x2) 
q{x) + q{xY 



q [x) 



(1 + 2x){q{x) + q{xY) - {q jx) + 2q{x)q {x)){x + x'^ 
q{x) + q{xY 



> 0. (93) 



derivative of p{y) is 
P'iv) = 



dp{y)/da 
dy/dx 

q'{x)- 



<0. 



' q{x) +q{xY 

The second-order derivative is given by ( |93l ), which is positive, as follows. 

Thus, p{y) is a convex function of y. ■ Proof of Proposition EJ 



5(//,X)} must be above the line segment connecting (si,S2) 
and (ii, ^2); and hence are not attainable by solving WSISMin. 
This is true for all /i's; hence if a boundary point is not at- 
tainable by solving WSISMin, it must be above a line segment 
connecting two particular points in S{fi,I) for some fi. With 
the above two lemmas, we are ready to prove Proposition |2] 



A 



(a) A straight line in X. (b) Tlie image of Q1Q2 in Ti- 

Fig. 12: Illustration of Lemma [T] Mapping a straight line in I to a 
convex curve in 7?.. 



Proof: First we define 
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(a) Inverse-SNR region 



(b) Rate region 



Fig. 13: Inverse-SNR region and corresponding rate and region. 



According to the above lemma, the line segment Q1Q2 in 
Fig. |12(a)| is mapped to a convex curve Q'iQ'2 in Fig .|12(b)| by 
U. In addition, it is easy to see that QiQ2 + li. i-> Q'iQ'2 ^ K, 
i.e., any point above Q1Q2 (for example, P in Fig. |12(a)| will 
be mapped to be a point below Q'iQ'2 (i.e., P in Fig. |12(b)| i. 

Lemma 2: Let a point (gi, 52) G bd(I) \ S{^,X). If qi = 
\ti + Xsi, where (^1,^2), (si, £2) G 5(/i,I) and (^1,^2) 7^ 
(si,S2), we have (72 > ^h + As2, i.e., the point (91,92) is 
above the line segment connecting (^1,^2) and (.si,S2)- 

Proof: We show this by contradiction. Suppose S{^,I) 
has more than one elements for a given /i, such that 
(ii,i2),(si,S2) e S{fi,I), and (^1,^2) 7^ (si,S2). Accord- 
ing to the definition of S{^,I) given by ( fTSl l. we have 
/^ti + /ii2 = /^Si + MS2 = fn^ where m is the minimum 
value of the weighted sum for a given fj, over all points in I. 
If qi = Ail + Asi and 92 < At2 + As2, we have 

IJ.qi+pq2 < fJ.{Xh+Xsi)+fi{Xt2 + Xs2) (94) 
= A(^ti +/2i2) + A(/^si +/is2) (95) 
= TO. (96) 

If fj,qi + Jlq2 < m, it contradicts that to is the minimum value 
of the weighted sum for the given fi; If fiqi + p,q2 = m, 
it contradicts that (<Zi,<Z2) is not in S{fi,I). Therefore, the 
lemma holds. ■ 

According to the above lemma, for a given /i, if S{fi, I) has 
more than one elements, the set of boundary points {((71,(72) ■ 
((71,(72) e hd{I)\S{fi,I),si < qi < ii,(.si,s2),(ti,t2) e 



A = {/i : 5(/i,T) has more than one elements}, 



(97) 



and let ^^ be the line segment (e.g., AB in Fig. |13(a)[ ) with 
two end points from S{i-i,I) for /i G A (e.g., points A and B 
in |13(a)| i. According to Lemma|2] the boundary points that are 
not attainable by solving WSISMin, denoted as bd(I) \ S{I) 
(here referring to curve AB in Fig. |13(a)| i, must be above 
l^'s, i.e., bd(I) \ S{I) C U^gaC'p + K); and it follows that 
U{hd{I)\S{I)) CU{\J^^^{lf, + K)). According to Lemma 
m ^iU^eAil^^ + K)) C U^eAml^.) - K), where Uil,.) is 
a convex curve (e.g., here U{lfj,) refers to the dashed convex 
curve A'B' in Fig. |13(b)| l. Let /^ be a line segment (i.e., the 
dot-dashed line segment A'B' in Fig. |13(b)| i that connects the 
two end points of the convex curve Z^(Z^). Due to the convexity 
of U{lf,), we have U{lf,) - K C f^ - K and hence U{hd{I) \ 
S{I)) C U^eA(^~M-K). Notice \J^^^^{l^+K) is sufficient for 
constructing O by convex hulling. Therefore, U {bd{I)\S {!)) 
or bd{I)\S{I) is not necessary for constructing O. Since B C 
bd(I), the set B\S{X) is also not necessary for constructing 
O. m 
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